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GAUSSIAN PROCESSES, KINEMATIC FORMULAE AND 
POINCARE'S LIMIT 
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We consider vector valued, unit variance Gaussian processes de- 
fined over stratified manifolds and the geometry of their excursion 
sets. In particular, we develop an explicit formula for the expectation 
of all the Lipschitz-Killing curvatures of these sets. Whereas our mo- 
tivation is primarily probabilistic, with statistical applications in the 
background, this formula has also an interpretation as a version of 
the classic kinematic fundamental formula of integral geometry. All 
of these aspects are developed in the paper. 

Particularly novel is the method of proof, which is based on a an 
approximation to the canonical Gaussian process on the n-sphere. 
The n — ► oo limit, which gives the final result, is handled via recent 
extensions of the classic Poincare limit theorem. 



1. Introduction. The central aim of this paper is to describe a new result 
in the theory of Gaussian related fields, along with some of its implications 
to both geometry, and to a lesser extent, to statistics. 

The basic object of interest is a M fc valued random field y defined on a 
parameter space M and its excursion sets 



where Dei' For most of the paper, we shall take M and D to be C 2 
stratified manifolds in M. N and respectively. 
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Stratified manifolds are basically sets that can be partitioned into the 
disjoint union of manifolds, so that we can write 

dimM 

(1.2) M= \_\ djM, 

3=0 

where each stratum, djM, < j < dim(M), is itself a disjoint union of a 
number of j'-dimensional manifolds. A typical 3-dimensional example is given 
by the saggy couch of Figure 1, in which case d%M is the interior of the couch; 
diM the collection of the six sides, some concave and some convex; d\M is 
made up of the 12 edges; and 8qM contains the 8 corner vertices. In Figure 
1, an excursion set might be the grey area where the stress y is greatest. 

Our aim is to study the global geometry of excursion sets, as measured 
through their Lipschitz-Killing curvatures, Cj {A{f, M; D)), j = 0, . . . , dim(M ) . 
In particular, since these curvatures are random variables, we shall be in- 
terested in computing their expectations. We shall define Lipschitz-Killing 
curvatures below. However, if you are unfamiliar with them, at this stage it 
suffices to know that £n(A) is a measure of the volume of A, Cn-i(A) a 
measure of its surface area, and Cq(A) its Euler characteristic, an important 
topological invariant. 

We cannot do this for all y. For a start, y will both have to be smooth 
enough for basic differential geometric techniques to be applicable. Thus, 
a basic requirement will be that y has, with probability one, C 2 sample 
paths. Furthermore, writing y = (yi, ■ ■ ■ ,yk), we shall assume that the yi 
are independent, identically distributed (hereafter i.i.d.) centered Gaussian 
processes of constant variance, which we take to be 1. The processes are, 
however, not assumed to be stationary. For such a y, we shall prove that 

dimAf— i r 



(1.3) E{d(A(y,M,D))} 



E 

3=0 



l+J 
j 



(27rri/ 2 £ t+j (M)M](D), 



where the combinatorial flag coefficients are defined below at (4.3) and the 
, described and defined in Sections 3 and 6.1, are certain (Gaussian) 
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Minkowski functionals that to a certain extent, play the role of Lipschitz- 
Killing curvatures in Gauss space. The Lipschitz-Killing curvatures on both 
sides of (1.3) are computed with respect to a specific Riemannian metric 
induced on M by the component processes yj. Note, however, that Cq{A) is 
the Euler-Poincare characteristic of A, and so independent of any Rieman- 
nian structure [cf. Theorem 4.1 for a formal statement of (1.3)]. 

1.1. What is new here? The result (1.3) has a long history. If M is a 
simple interval [0,T], y is real valued and stationary, and D = [u, oo), then 
(1.3) is essentially the famous Rice formula, which gives the mean number 
of upcrossings of the level u by /, and dates back to 1939 [15] and 1945 [16]. 
Since then, there have been tens, if not hundreds, of papers extending the 
original Rice formula in many ways, with the developments up until 1980 
summarized in [1]. More recently, there was a series of papers by Worsley 
(e.g., [24, 25, 26, 28]) that were important precursors to the general theory 
presented in this paper. However, the first precursor to (1.3), at the level of 
processes over manifolds with C 2 boundaries, appeared only in 2002 in [21], 
where we considered only the first Lipschitz-Killing curvature £o(A u (f, M)) 
and then only for real valued y. In [20], one of us (JET) extended this to 
vector valued y, which allowed for the derivation of the far more general, 
and far more elegant, (1.3) for the first Lipschitz-Killing curvature. 

What is new here then is the extension to parameter spaces as general as 
stratified manifolds, and the extension to all Lipschitz-Killing curvatures. 
Both of these are important for applications. However, perhaps more impor- 
tant, and certainly more novel than either of these, is the method of proof. 
The proofs in the current paper are new, and far more geometric than the 
earlier ones. In particular, the proof in [20] progressed primarily by evalu- 
ating both sides of (1.3) and then showing that they were equivalent. The 
current proof starts on the left-hand side and, eventually, yields the right- 
hand side. The geometric nature of the current proof also explains why the 
two sides should be equal. 

1.2. Statistical implications. The general structure of (1.3) has signifi- 
cant implications for a class of problems out of the purely Gaussian sce- 
nario. Taking F :M. k —> M to be piecewise C 2 , along with appropriate side 
conditions, and defining a (now non-Gaussian) process 

(1-4) f(t)=F(y(t))=F( yi (t),...,y k (t)), 
with y Gaussian as above, it follows immediately from (1.3) that 
E{£i(A(f,M, Koo)))} 

(1.5) 

dim M—i r . 

= £ 1 + 3 (2^/ 2 £ i+j (M)M](F- 1 [u,+oo)). 

3=0 L J J 
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Non-Gaussian processes of the form (1.4) appear naturally in a wide vari- 
ety of statistical applications of smooth random fields (e.g., [2, 3, 4, 19, 24, 
25, 26] with an excellent introductory review in [27]). 

An additional and extremely important application of (1.3) lies in the so 
called "Euler characteristic heuristic" that for a wide range of random fields 
/, 

p{sup/(i)>u)-E{£ (4(/,M,[u,oo)))} < error(u), 

where error (u) of a smaller order than both of the other terms as u — > oo. In 
the Gaussian case, this heuristic is now a well-established theorem, and the 
error term is known to be of order exp(— u 2 (l + rj) /2) (for an identifiable rj > 
0) while both the probability and expectation are of order exp(— u 2 /2) [22]. 
The ability to compute the expectation therefore provides useful, explicit 
approximations for the excursion probability. 

1.3. Geometry. One of the basic results of integral geometry is the so- 
called kinematic fundamental formula (henceforth KFF), which in its sim- 
plest form, states that for nice subsets M\ and M 2 of W 1 , 



(1.6) 



C i (M 1 n g n M 2 ) dv n (g n ) 



3=0 







n 


i 




J. 



C i+j {M x )C n ^{M, 



2 ■ 



Here, G n is the isometry group of W 1 with Haar measure v n normalized so 
that for any x £ W 1 and any Borel A C M. n , v n ({g n £ G n : g n x £ A}) = H n (A) , 
where 7i n is n-dimensional Hausdorff measure. (See [12, 18] for Mj elements 
of the convex ring or similar, and [6] for more esoteric Mj closer to the spirit 
of this paper.) 

Now reconsider (1.3). Taking (f2,.F,P) as the probability space on which 
y lives, (1.3) can be rewritten as 



(1.7) 



d(M D{y(Lo)y 1 D)dF{u J ) 

(2Tr)-i/ 2 £ i+j (M)M](D). 



dim M—i 
3=0 



l + J 

3 



Written this way, it is clear on comparing (1.6) and (1.7) that our main 
result can now be interpreted as a KFF over Gaussian function space, rather 
than over the isometry group on Euclidean space. We find this interpretation 
novel and intriguing bridging as it does between a probabilistic problem and 
a geometric answer of classic form. 
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1.4. More on stratified manifolds. Although (1.2) gives the basic struc- 
ture of a stratified manifold, certain assumptions need to be made for the 
results of this paper to hold. In particular, we need to assume that each 
connected component in each stratum is a C 2 manifold. More importantly, 
we need to assume that both M and D are Whitney stratified spaces (see 
[10, 14]) which ensures that the various strata are "glued together" in a 
smooth way. There are further technical assumptions of "tameness" and 
"regularity" that can be found in Chapter 15 of [3]. One of these is local 
convexity; a demand much weaker than convexity (the saggy couch is locally 
convex) , which while not necessary, is something that we shall assume in this 
paper so as to make many formulae somewhat more manageable. Handling 
all these technicalities is not only tiresome, consuming both space and time, 
but often also quite difficult. Nevertheless, we shall completely avoid them 
since our main aim in the current paper is primarily didactic. We want to 
present the main result in a setting that gives the flavor of the general re- 
sult, introduces the new style of proof and yet does so in a way that is both 
self-contained and comparatively easy to follow. If you want all the details, 
then Chapter 15 of [3] is the place to turn. If you want rigor without leaving 
this paper, then you can assume that M is a manifold without boundary 
(necessarily locally convex) such as a sphere or a simple cube, which has 
all the edges and corners but none of the curvature of Figure 1. In neither 
of these cases are additional assumptions needed, and for both of them the 
results and proofs of the paper are still new. 

1.5. A plan of action. Our first task, in the following two sections, will 
be to define the Lipschitz-Killing curvatures Cj and the Gaussian Minkowski 
functionals Mj. Once this is done, then at very least all the terms in (1.3) 
will be clearly defined, and we shall be able to state the formal result as 
Theorem 4.1. Thus, the reader who is not at all interested in proofs can 
leave us at that point. 

We shall then slowly start on the proof, proving in detail only a special 
case of the Theorem 4.1 for the so-called isotropic process on the unit sphere 
in M. N . This is the most interesting part of the paper, and contains not only 
new results, but more importantly, a conceptually new way of looking at 
things. Despite limiting the proof to one special case, we shall explain in 
Section 8 why this case is particularly important, and also why all other 
cases can, in principle (but with quite a lot of additional work) be shown to 
follow from it. 

However, even for this specific choice of process, the computations are far 
from trivial and so we approach them in an indirect fashion. In particular, 
we shall start by looking at non-Gaussian processes with finite expansions 
and with coefficients coming from random variables distributed uniformly 
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over high dimensional spheres. This will enable us to take expectations using 
an appropriate version of the KFF on spheres. 

The passage from this scenario to the Gaussian one will be via a limit 
theorem for projections of uniform random variables on the n-sphere as 
n — > oo, historically associated with the name of Poincare, although we shall 
need a slight extension of some more recent versions due to Diaconis and 
Freedman [7] and their generalizations to matrices in [8]. 

2. Lipschitz Killing curvatures. Perhaps the easiest way to meet Lipschitz- 
Killing curvatures is via Weyl's tube formula, which in our setting states that 
for quite wide classes of sets M in M. 1 , and for sufficiently small p > 0, there 
exist numbers Cj(M) for which 



where Tube(M, p) = {t E M : inf s& M \ s — t\ < p}, and Lo m is the volume of the 
unit ball in M m [cf. (4.3)] and, in a notation anticipating something more 
general, Hi is /-dimensional Lebesgue measure. 

Weyl's formula actually defines the Cj(M) for this Euclidean case for all 
< j < dim(Af), and we define Cj{M) = for all j > dim(M). However, 
they can also be computed directly, and have a natural extension to the 
situation in which M is a Riemannian manifold, a situation which we shall 
soon need, even though our parameter spaces M are basically Euclidean. In 
the Riemannian case, we shall need to assume that M is embedded in an 
ambient manifold M. We shall denote the Riemannian metrics by g and g, 
respectively, the curvature tensors of M and M by R and R, and denote the 
scalar second fundamental forms of the strata djM as they sit in M by 5. 
We also have Riemannian (Hausdorff) volumes on M and djM, which we 
continue to denote by TCn and 7ij, respectively. 

In this more general setting, Lipschitz-Killing curvatures of M are defined 



dim M 



(2.1) 



H l (Tube(M,p))= Y, ^i-AW, 



by 



dim M 



Ci{M) 



J2 (2vr)-^- i )/ 2 



j=i 



LCJ-0/2J 




2m) 



x E 



(2.2) 



m=0 




j—i—2m 



) 



x l NtM Hi-j-i(drj)Hj(dt), 
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where TtdjM is the tangent space to djM at t, TtdjM is the normal 
space, N t M is the normal cone to M at t, Tr Tt9jM is the trace in T t djM, 
and S(TtdjM ) is the unit sphere in TtdjM . [We shall generally write 
S\(V) and -Ba(^) for the sphere and ball of radius A in a vector space V, 
with S = Si and B = B\.\ The constants C(m, i) are given by 



and s n = r(^)/27r n//2 is the surface measure of 5(R n ). For i > N, we set 
£j = 0. We shall always assume that Lipschitz-Killing curvatures are finite. 

If you decided earlier to take the route of avoiding the generality of strat- 
ified manifolds, and you want to think of M as being a smooth manifold 
without boundary, and you want to work with the standard Euclidean met- 
ric, then note that there is only one element in the outer sum of (2.2), which 
is the case j = N. This makes things much easier. Alternatively, if you want 
to think of M as the simple, Euclidean, A r -dimensional rectangle ni=i [0> ^i] ; 
then both it! and S are identically zero, and it is a simple exercise to see that 
£j(M) = J £2Ti 1 - ■ -Tij, the sum being taken over the ( .) distinct choices of 
subscripts. 

For further details, such as a proof of the fact that Cj are independent of 
the stratification; see, for example, [3, 5, 11, 23]. 

The fact that Cq(M) is equivalent to the Euler-Poincare characteristic 
of M, and so independent of any Riemannian structure, is the celebrated 
Chern-Gauss-Bonnet theorem. The remaining Lipschitz-Killing curvatures 
also appear under a variety of other names, such as quermassintegrales, 
Minkowski, Dehn and Steiner functionals, and intrinsic volumes, although 
in many of these cases the ordering and normalizations are different from 
ours. 

We now almost have enough to decipher the meaning of the Lipschitz- 
Killing curvatures appearing in the main result (1.3). What remains is to 
define the Riemannian metric under which they will be computed for most 
of this paper. 

Recall that Riemannian metrics are defined on pairs of vectors in the 
tangent spaces Tf. Thus, taking Xt,Yt€ TjM, letting yi : M — > R to be one of 
the components of our random function y, and writing X^y for the derivative 
of y at t and in direction Xt, we can define the (Riemannian) metric induced 
by y as 



where C(s,t) = E{j/j(s)yj(f)} is the common covariance function of the yi. 
(We implicitly assume that C is such that g is truly a Riemannian metric, 
in that gt is nondegenerate for each t £ M.) 




m + i > 



m = 0, 



(2.3) g(X,Y) = g t (X t ,Y t )^E{(X m ) • (Y m )} = X s Y t C(s, t)\ s=t 
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Note from (2.3) that all the tools of Riemannian manifolds, connections, 
curvatures, etc., can be expressed in terms of covariances. In particular, 
it turns out that all of these tools also have interpretations in terms of 
conditional means and variances. See [3] for details. 

One case, however, is worth noting now. If y is locally isotropic, with unit 
variance derivatives in all directions, in the sense that 

K{X m (t)Y m (t)} = X s Y t C( S ,t)\ s=t = {X t ,Y t ), 

where the inner product is the usual one, then the Riemannian metric is the 
usual Euclidean one. 

3. Gaussian Minkowski functionals. Having determined the meaning of 
the Cj(M) of the main result (1.3), we now turn to the A4j(D). Again, a 
tube formula approach is the easiest to take, at least initially. For this, let 
7 = Ik be Gauss measure on M. k , so that for Borel icR* 1 , 

lk (A) = (27r)- k / 2 [ e-\ x \ 2 ' 2 dx. 
J A 

Once again, take Del* to be a stratified manifold, satisfying the same 
conditions we set out for M. We also need an an additional integr ability 
condition, a weak form of which is given in Chapter 15 of [1], but which we 
shall bypass by assuming that all (Euclidean) curvature tensors and second 
fundamental forms of D and its strata are uniformly bounded. 

Then there exist numbers .Mj(M) = A4j fc (M), j > 1, called the Gaussian 
Minkowski functionals so that, for p small enough, 

OO A 

(3.1) 7*(Tube(Ap)) =lk{D) + ]T ^Mf(D). 

3=1 J ' 

This gives an implicit definition of the A4j , akin to that which Weyl's tube 
formula (2.1) gave for the Cj. We shall give a more constructive definition, 
in the spirit of (2.2), in Section 6.1. 

In many cases, the A4j are quite easy to compute. For example, if M = 
[it, oo) C M, then since Tube([u, oo), p) = [u — p,oo), comparing a Taylor se- 
ries expansion of 7(Tube(M, p)) with (3.1) easily gives, for j > 1, 

Mf([u,oo)) = {2i:y 1 ' 2 H^ 1 {u)e- u2 l 2 , 

where H n is the nth. Hermite polynomial. Additional examples, also acces- 
sible from simple calculus, can be found in [3, 20]. For a proof of (3.1), see 
either [20] where, to the best of our knowledge, it first appeared for manifolds 
with C 2 boundary, or Chapter 10 of [3], which treats stratified manifolds. 
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4. The main result. We now have all that is needed for understanding 
the result (1.3). Indeed, with just one more condition, we can state it prop- 
erly. 

Writing for the moment y to represent one of the i.i.d. real valued com- 
ponents of our W k valued random process y, assume that it is Gaussian, 
centered, of constant variance one and possessing almost surely C 2 sample 
paths. Assume furthermore that for each t € M, the joint distributions of 
{y,dy/dti,d 2 y /dtj dtk)ij,k=i,...,N at t are nondegenerate, and for some finite 
K, and all s,teM, 

(4.1) max\C ij (t,t) + C ij (s,s)-2C ij {s,t)\ < K\]n\t - s\\~ (1+a) , 

where Cjj = d 4 C/d 2 tid 2 tj is the covariance function of d 2 y /dtidtj. 

Theorem 4.1. Let M C R N and Del* be stratified manifolds satisfy- 
ing the conditions described in Sections 2 and 3. Let y = (yi, . . . ,yk) '■ M — > 
M fc be a vector valued random process, the components of which are indepen- 
dent, identically distributed, real valued, Gaussian processes satisfying the 
above conditions. Then 

dim M—i r . .-, 

(4.2) E{C l (Mny- 1 (D))}= £ l + J (2^/ 2 C i+j (M)M](D), 

3=0 1 J J 

where the Cj, j = 0, . . . ,N are the Lipschitz-Killing measures on M with 
respect to the metric induced by the yi, the A^J are the Gaussian Minkowski 
functionals on R fe and the combinatorial terms are given by 



[n\\ = n\uj n , lo ti 



(43) [k\ [k]\[n-k]lu; n cu n _ k > L '" J - ^ r(f + l)' 

Lo n being the volume of the unit ball in W l . 

Now that we know what all the terms above mean, it is worthwhile to look 
again at the structure of (1.3) / (4.2). Note that as far as the right-hand side of 
the equation goes, there is an elegant factoring of parameters. The covariance 
structure of y appears only in the Cj(M), as does the geometry of the 
parameter space, while the M'J(D) are computed without any reference to 
either of these. This is particularly useful for applications such as (1.5), since 
it implies that when working with / of the form / = F(y) the computations 
related to M and the distributional properties of y separate out completely 
from the properties of F. 

We can now turn to the proof of Theorem 4.1, which will take up the rest 
of the paper. In fact, we shall give a full proof for only one, very specific 
process, on a very specific parameter space. This special process is defined 
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Fig. 2. The pre-Gaussian process from a subset of a great circle to R 2 . 

on unit sphere S(W) and is known as the canonical isotropic process. Its 
components are independent, centered Gaussian processes with common co- 
variance function 

(4.4) E{ yi (s) yi (t)} = (s,t), 

where (•,•) is the usual Euclidean inner product. You can jump forward 
to Section 8 to see why this process is so special, and why all other cases 
effectively follow from it. Now, however, we shall begin proving Theorem 4.1 
for iV-dimensional M C S(M. 1 ), l> N, and for the canonical process. 

In fact, even the canonical isotropic process will be a little too hard for us 
to work with directly, so we shall approximate it by something even simpler. 
The approximations are given by a sequence, {y^ n '}n>h of smooth M fc -valued 
processes on S(M. 1 ). To define them, for each n > Z, we first embed S'(M') in 
S"(IR n ) in the natural way, by setting 

S(R l ) = {t = (t 1 ,...,t n )e S(R n ) : t l+1 = -.. = t n = 0}. 

Taking the rotation group 0(n), equipped with its normalized Haar measure 
u n , as an underlying probability space, the nth process j/ n ) is defined by 

(4-5) y {n \^9n) = ^^i, n , h {Vngnt), 

where t £ S(R l ), g n £ 0{n) and n^ n k is the projection from S^(R n ) to 
]R fc given by 

(4.6) vr A>nifc (xi,...,x n ) = (x!,...,x k ). 

Figure 2 gives an example of this when M is a subset of a iV = 1 dimensional 
great circle in the 2-sphere (1 = 2), taking values in the plane (k = 2). 

Now recall the old and quite elementary result known as the Poincare limit 
theorem, or Poincare's limit (although whether it really is due to Poincare 
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is not clear [7]). In our current setting, a recent more powerful version due 
to Diaconis, Eaton and Lauritzen [8] implies that the finite dimensional dis- 
tributions of ?/ n ) converge, in total variation norm, to those of the canonical 
process on R'. More importantly, for functionals F of these processes for 
which E{|F(y)|} < oo, we have 

(4.7) lim E{F(y^)} = E{F(y)}. 

Thus, if we can compute mean Lipschitz-Killing curvatures for the approx- 
imating processes, we shall be in good shape obtaining the result we are 
seeking by a limiting argument. The key to the first step lies in a version of 
the KFF on spheres. 

5. The KFF on S\(R n ) and the beginnings of a proof. Although we met 
the KFF for Euclidean spaces back in the Introduction [cf. (1-6)], we need 
a version for the sphere 5(R n ), good references for which include [9, 13, 17, 
18]. To formulate it, we need to extend the Lipschitz-Killing curvatures of 
Section 2 to a one parameter family of measures Cf, for k > and < i < 
N = dim(M), by replacing the Riemannian metric R in (2.2) by R + (k/2)I 2 , 
where / is the identity tensor. As for the basic Lipschitz-Killing curvatures, 
for i > N, we define Cf = 0. 

We also define a one parameter family of curvature integrals Cf(M) = 
Cf (M, M) . For lack of better terminology, we call these new curvature mea- 
sures extended Lipschitz-Killing curvatures. It is trivial that C® = C{. It 
takes some algebra (cf. [3]), but it is not too hard to show that there are 
simple equivalences between the Cj and CJ, given by 

^ ^_y(-^) n (i + 2n)\ ~ K n (i + 2n)\ 

(5.1] L,,-^ L i+2n , Ci-2_, ,., L i+2n . 

n=0 v ' n=0 v ' 

[Note that the summations here are actually finite, since all curvatures are 
zero for i > dim(M).] 

It turns out that in dealing with subsets of «S , A(R n ) the Cj are more 
natural to deal with than are the original Lipschitz-Killing curvatures. For 
example, there is a elegant version of Weyl's tube formula involving them, 
and, more importantly for us, there is a nice KFF. 

Let Gn t \ denote the group of isometries (i.e., rotations) on S\(R n ), with 
Haar measure v n \ normalized so that for any x £ S\(M. n ) and every Borel 
A C S x (R n ), v n ,\{{gn 6 G n y.g n x £ A}) = H n -i(A), where U n -\ is surface 
(Hausdorff) measure. The KFF on S\(R n ) then reads as follows, where M\ 
and M<i are stratified manifolds in S\(W n ) satisfying the conditions that we 
have been assuming: 

I Cf 2 (M 1 ng n M 2 )du n>x (g n ) 
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(5.2) 



E 

ji— 1— i 
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-1 

••A 







n — 1 


i 




3 



-£^.(M 1 )£^ W (M 2 ) 



Now let us apply this KFF to the processes y^ n \ The key result, from 
which everything else follows, is the following, under the usual conditions 
on M and D. 

Lemma 5.1. Let y^ be the model process (4.5) on M C S(R l ), with 
n>l. Then for D cR k , 



dim M—i 

(5-3) = £ 

3=0 





n — 1 


1 


+ 




. 3 „ 




. 3 . 



1 / n r\ ■? \/n,n,k ' 



C +i (M) 



SnU^" 1 )/ 2 



dim A4"— i 



^ S i+ i x ^n-l-jl^fc.n.fc 1 ^ 



[It is important to understand the meaning of vr^i- n k D above, and in all 

that follows. The problem is that for all t £ S ^(R n ), vr ^ >n>k (t) G B ^(M n ), 
which may or may not cover D. Thus, since 

it follows that tt~}~ , D may be only the inverse image of a subset of 
D.] 

Proof of Lemma 5.1. Since 7r"i ,D is a C 2 domain in S(IR n ), it 
follows from the construction of y^ n > that 
E{Cl(Mn(y ( - n) )~ 1 D)} 



0(n) 



4(M n n-V2 5 -i (7r ^_ n kD)) dfln{gn] 



O(n) 

n^/ 2 / £r i ( v ^Mn 5 - 1 (7r^ nfeJ D))^ n ( 5n ) 



O(n) 



(n-l+i)/2 



G. 



, y/n 



C> (^MC\g n {^ nk D))dv n ^{g n ), 
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where the second-to-last line follows from the scaling properties of Lipschitz- 
Killing curvatures and the last is really no more than a notational change, 
using the definition of y n ,\- 

However, applying the KFF (5.2) to the last line above, we immediately 
have that it is equal to 



dim M—i 

E 

3=0 



i/2 



~i + f 




n — 1 


i 




3 



n (i+i)/2 



= E 

which proves the lemma. □ 



n^ 2 


n — 1 


-l 


+ 




3 




3 



^-l-3^n,k D ) 



s n n 



(n-l)/2 



Suppose we send n — ► oo in (5.3), which by Poincare's limit is effectively 
equivalent to replacing the model process with a M fc valued canonical 
Gaussian y. Then in order for E{Cj(M C\y~ 1 D)} to be finite for the limit 
process y, we would like to have the following limits existing for each j < oo: 



(5.4) 



lim n 

n— >oo 



3/2 



n — 1 

i 



- 1 ^X-i^njP) 



Snn (n-l)/2 

A Stirling's formula computation shows that if the limit here exists, then 

-l 

(5.5) p j (D) = (2ir)-i/ 2 \j]\ lim ' 



Snn (n-l)/2 

Sending n — > oo in Lemma 5.1 and applying Poincare's limit (4.7), we see 
that if EUCKMny^D^} < oo, then 



K{Cj(Mny- 1 D)}= lim E{C}(M n {y {n) )~ l D)} 



(5.6) 



dim M- 

E 

j=0 



i+j 
i 



This is starting to take the form of (1.3) and (4.2). The combinatorial flag 
coefficients are in place, but both sides of the equation are based on the 
curvatures rather than the and we have yet to identify the functions 

Pj. Note the important fact, however, that on the right-hand side of the 
equation we have already managed a split into product form, each factor of 
which depends on the underlying manifold M or the set D, but not both. 

Moving to the desired curvatures, under the assumption that the limits 
(5.4) exist is the next step, and the result is summarized in the following 
lemma, again with the usual assumptions on M and D. 
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Lemma 5.2. Let M c S(R l ) and assume that for < j < dim(Af), the 
Pj(D) of (5.5) are well defined, finite and 

(5.7) E{\d(M ny~ 1 L>)|}<oo 

for y the M. k valued canonical isotropic Gaussian process on S(R l ). Then 

dimAf— i r 



(5.8) 
where 



E{£ ! (Mn^ 1 D)}= J2 

1=0 



i + 1 
I 



LC?-i)/2j 



Pj(D) 



(7-1)1 E 



^ (47r)'Z!(j-l-20! 



£+,(M)fl(£>), 



Pj-2l(D), 



i>i, 

j = 0. 



Proof. As usual, set iV = dim(M). Combining (5.1) and (5.6), we have 
E^ifMnf 1 !))} 



L(JV-i)/2J 

E 

m=0 



1 (, + 2 *"> ! -4 J ,(Mn ! ,-. D )' 



(47r) m m! i! 



L(A V /2J 1 (i + 2m)! JV -^; 2m 
^ (4vr) m m! i! ^ 



m=0 v / ' jf=0 

1 (i + 2m)\ 



i + 2m + j 
j 



C 



l 

i+2m+j 



(M)pj(D) 



L(jv-0/2J 

E 

m=0 



(47r) m m 



i! 



N-i-2m r 

x E 

5=0 



z + 2m + j 
j 



Pj(D) 

l(N-i-2m-j)/2\ 



E 

1=0 



(-I) 1 (i + 2m + j + 2l)\ 
(4tt) z Z! (i + 2m + j)! 



The passage from the above to (5.8) is now one of algebra, and the details 
of the three page calculation can be found on pages 408-410 of [3]. Since 
they add no insights, and since our aim in this paper is mainly to get the 
basic ideas across, we leave them to you to check. □ 

What remains now to do is to check that the limits (5.5) are indeed well 
defined, and to evaluate them. Before we can do this, we need to make a 
small excursion into differential geometry. 
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6. A geometric interlude. We need to better understand Gaussian 
Minkowski functionals and, in particular, their relationship to Lipschitz- 
Killing curvatures. In the following section, we shall use this information to 
evaluate the limits (5.5), and via them the all important functions Pj(D) of 
(5.9). 

6.1. Gaussian Minkowski functionals. We now give an explicit constru- 
tion of Gaussian Minkowski functionals, rather than relying on their implicit 
definition via the tube formula of (3.1). First, however, we need to somewhat 
extend the notion of Lipschitz-Killing curvatures. 

We shall continue to assume that M is a iV-dimensional stratified manifold 
in R J , locally convex and with uniformly bounded curvature. Take Borel A C 
R l and B C S(M. 1 ), retain the notation of (2.2) and define for < i < I - 1, 
a family of generalized Lipschitz-Killing curvature measures supported on 
M x S(R') by 



Ci{M,A x B 

limM 



dim A/ 



j=i 

LCJ-0/2J 
(6.1) x £ 



m=0 



-l) m C(l-j,j -i-2m) 
m\(j — i — 2m)! 



f f Tr T t d 3 M {R m Si~y 

JdiMnAJs(T ± d i M)nB 1 3 



IdjMnA JS(T ± d j M)nB 

x l NtM Hi-j-i{dui-j)Hj(dt). 

For i = I, we define /^j only on sets of the form A x S(R. 1 ) by setting Ci(A x 
S'(R')) = Hi(A). For Borel / : R l x 5(R Z ) -» R let £,(M, /) denote the integral 
of / with respect to Hi. 

A change of numbering and normalization now defines the generalized 
Minkowski curvature measures as 

(6.2) Mj(M, A x B) = {j\uj)Ci^{M,A x B). 

With these definitions, we can now give a direct definition of the Gaussian 
Minkowski functionals appearing in the tube formula (3.1) by setting 



(6.3) M](M) 4 (2^/ 2 £ f J m 1 )^ m+1 (M,^ 1 „ m ((t,r ? ))e-l t l 2 /2 



m=0 



Here, is the nth Hermite polynomial and (t,rj) is the standard inner 
product between i£M and 77 a vector in S(M. 1 ). 
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6.2. Warped products of Riemannian manifolds. Note first that the set, 
7T ^ nk D 6 S^(M. n ), so crucial to the arguments of Section 5 and crucial 
for all that will follow, can be stratified topologically into the disjoint union 
of the (n — k — 1 + j)-dimensional strata 

(6.4) D n ^ l+j ^D 3 xS(M n - k ). 

Furthermore, each such D n _ k -i+j can be written, again topologically as a 
disjoint union 

(6.5) D n _ k _ 1+j ~ (Dj n S^(]R fc )) U (Dj n (B Rk {0,V^))° x 5(M"- fc )). 

(See the example in Figure 3.) 

Since each Dj n S' v ^;(IR' c ) is a stratified subset of S^^R* 1 ), its Lipschitz- 
Killing curvatures can be computed using the tools we have developed so far. 
The second set in the union is, however, somewhat more complex. Although 
we have written it as a product set, this is only correct up to topological 
equivalence. We, however, need precise Lipschitz-Killing curvatures, which 
with the exception of Co, are not topological invariants. To handle this, we 
need we need to break the Riemannian structure of products into a product 
of structures. Second, we need to treat each such product, in which each copy 
of S(M. n ~ k ) is likely to have a different radius, as a subset of a Riemannian 
warped product. 

Recall that the Riemannian warped product of two Riemannian manifolds 
(Mi,gi) and (A/2,52) with a smooth warp function a 2 : Mi — ► [0, +00) is the 



/ 





u \ 


n V2 









Fig. 3. A particularly simple example of (6.5) when N — k — 1, n — 3, and D = [u, 00) = 
Dq U Di = {u} U (u, 00). In this case, 7r^i g 6 Sy^R 3 ) is the spherical cap shown. For 

j = 0, Z?3_i_i_|-o = Di is the set made up of the left-hand boundary of the cap along with 
its rightmost point. For j — 1, D2 is the rest of the cap. 
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Riemannian manifold 

(6.6) (M 1 ,M 2 ,a) = (Mi x M 2 , 9l + a 2 g 2 ). 
Usually, as in our case, M 2 is a sphere. As an example, consider 

M a = (B Rk (0,V^))° xS(M n - k ), 
where the Riemannian metric on the open ball is given by 

(6.7) g a = g Rk + Va Va 

for <7 2 (t) = n — ||t||gfc, and the Riemannian metric on S(M. n ~~ k ) is the canoni- 
cal one inherited from M. n ~ k warped by o~ 2 (t). The importance of this exam- 
ple to us is that each stratum D n ^k~i+j is isometrically embedded in such 
a warped product. Using this embedding, we shall be able to compute 

,1/n 



C n-l-i(*Ji, n ,k D > D n-k-l+j) 



the contribution of these strata to the Lipschitz-Killing curvatures of ir^ n fe D. 

6.3. Connections and second fundamental forms. The first step to com- 
puting these contributions is to determine the form of the Levi-Civita con- 
nection V 7 of M a , as this is needed to in order to compute the second 
fundamental form of Dj x S(W l ~ k ) in M a . Therefore, consider a general 
warped product (Mi,M 2 ,a) and denote the Levi-Civita connection on each 
Mj by V J . Use E, or Ej, to denote vector fields on M±, identified with their 
natural extensions on M± x M 2 . Similarly, F or Fj, denote vector fields on 
M 2 extended to M\ x M 2 . Then from the definition of Levi-Civita connec- 
tions, the product structure of M and Koszul's formula is not too hard to 
check that 

V El E 2 = V El E 2 , V Fl F 2 = V 2 Fl F 2 - a 2 g 2 {F x , F 2 )V a 2 , 

(6.8) 

W E F = V F E = E{\oga)F. 

With the Levi-Civita connection defined, the next step lies in determining 
the second fundamental form over the sets Dj x S^M* 1 ) as they sit in M a , 
as well as traces of their powers. For this, we need to describe the normal 
spaces T^^Mq for (t,rj) £ Dj x S(M. n ~ k ). A simple argument shows that at 
these points, 

(6.9) (T^M*) 1 = (T t Dj © T ?) 5(M n - fc )) ± ~ T t Df, 

where TtDj- is the orthogonal (with respect to g a ) complement of TtDj in 
T t B Rk (0,^/n). 

We can therefore state: 
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Lemma 6.1. Retaining the above notation, for < I < n — 1, take 
(t, V ) € Dj x S(R n ~ k ), Vk _j € (T( t>r) )Dj x 5'(M ri_fc ))" L . 

Then 

r=0 

where S is the second fundamental form of Dj x S'(M n_fc ) in M a and S a is 
the second fundamental form of Dj in (B-^k(0,y/ri),g a ). 

Proof. Fix an orthonormal (under g) basis (E±, . . . , E k ,Fi, . . . , F n _ k _i) 
of T( t ^M a such that (E\, . . . ,Ej) forms an orthonormal basis of TtDj. Ob- 
servation (6.9) implies that any v k -j can be expressed as 

n— fc— 1 

Vk—j — ^ a r F r 

r=l 

for some constants a r . Applying (6.8) and the Weingarten equation, we find 

Sv k _ J {E r ,E s ) = —g(Vz r v k _j,E s ) = —g a (V^ r u k _j,E s ) = S a>Uk _. (E r , E s ), 
S Uk ^{F r ,F s ) = -g{V° Fr v k _j,F s ) = -y k _j{\ogc7 t )g{F ri F s ) 

= -u k _j(loga t )5 rs , 
S Vk _ j (E r ,F s )=0. 

Therefore, for each v k -j, the matrix of the shape operator S Uk _ j in our 
chosen orthonormal basis is block diagonal with one block, of size j, given 
by {S a ,v k _ • (E r , E s )}i< rtS <j and the other, of size n — k — 1, being given 
by — v k -j(logat)I(n-k-i)x(n-k-i)- Therefore, applying basic combinatorial 
properties of the trace operator we have that for I < n — k — 1, 

l^(Sl k J = £ ( n - k _- r l ) (-l)^(^(log^y-iTr(^.), 

r=0 

which completes the proof. □ 

7. Back to the main proof for the special case. With the geometry be- 
hind us, we now turn to the asymptotics required for computing the limits 
(5.5), for which the following lemma is a crucial step. The final move to the 
Gaussian Minkowski functionals appearing in Theorem 4.1 above will then 
involve no more than some careful asymptotics. We start with a lemma, 
assuming throughout that the conditions of Theorem 4.1 hold. 
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Lemma 7.1. Let D n _fe_i +J be as defined at (6.4) and i>k—j > 0. Then 



i+j — k 



=0 

x£ J _ i ( J D CT , C r" +fc - 2i - 2 ^ 2Z - 1 (27r)- fc / 2 / li+J _ fc _ i l Dj ), 

where D a = D n -B R fc (0, y^n) is i/ie regular stratified manifold obtained from 
the intersection of the embedding of D in M. k and the open ball of radius s/n, 
endowed with the metric g a given by (6.7) and hi(t,u) = (u,t)L k . 
Furthermore, with Lp k denoting the k- dimensional Gaussian density, 

£™ s n{ n-l)/2 C n-l-i^~An,k D i D n -k-l +J ) 

(7-2) 

= H ][k+ [i\\ ~ (fc+1- j- l) ^J-l( D ^k(t)h i+j - k -il Dj ). 

Proof. Note first that we can assume that n is large enough so that 
(7.3) Dn-k-1+3 n S^(R k ) = 0. 

Suppose this were not the case. Then we can further stratify Dj into two 
parts, one of which has no intersection with S^(M. k ) and one of which is 

contained in S^(M. k ). It is only the second part that concerns us. How- 
ever, for this part, we also have that 7r~i- , Z?,- = D,-, and so 7r~J- ,Dj is 

r 1 vn,n,fe J ■>' y/n,n,k ■> 

a 7-dimensional stratum of 7r _ i , D. However, such strata only contribute 

J • v rn,n,k J 

to Lipschitz-Killing curvatures of order to j. Since we are interested in 
(asymptotically) high-order Lipschitz-Killing curvatures, we can therefore 
forget these components and assume that (7.3) holds. 

Thus, from the definition of the Cf, and for alH > k — j > 0, we now have 

£i-l-i(\4n )fc A Dn-k-l+j) 

C(k-j,i + j-k) 



(2 7 r)( i +J- fc )/ 2 (i+j- jfe)! 

xS(m. n - k )JS(T (Uv) D :j xS(R n - k )- L ) 3 

x ± NtM H n -i-j+k(dt,dr]), 



where 



7U-i-j+k(.dt,dn) = <y r r x ' k H n -i-kidn)Hj{dt) 
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is the Hausdorff measure that Dj x S(M. n ~ k ) inherits from D a , the warped 
product of (D a ,g a ) and »S(R fc ) with its usual metric and warp function a 2 
as in (6.7). Then by Lemma 6.1 and (6.9), 

£ i-l-i( 7r VH,n,fc- D ' Ai-fc-l+j) 

= C(k + k){2Tr)-( l+ i~ k V 2 

i+j-k 



•sr-^ I n — 1 — k 

™- k )Js{T t Df) ^ \i+j-k-lj ut 



^n—l—k 
a t 



i+j-k-lf,. rr.Wi+j-k-l 



x (_ 1 ).+ J -*-l (l/fc _ i(logo . t) ) 

Xj^{S l ^ k _.)n k ^{du k ^) 

xHn-i-kidvWjidt). 
Equation (7.1) now follows from the fact that 

C(k-j,i+j - fc)(27r)-(^- fc )/ 2 _ sj+^j 
C{k-j,l){^)- l/2 Si ' 

followed by integrating over S(M. n ~ k ) and noting that 

v k -j (log a t ) = J ~2 • 

a t 

As for the second conclusion of the lemma, (7.2), note that 

S n -k ( n — k — 1 \ n+k~2i-2j+2l-l 
n^o Snn (n-l)/2 y i+ j-k-l) at 

= Hm f. " T ^ 7 1 „(n+fc-2i-2i-2i-l)/2 

Snn (n-l)/2 y l+J -k-lj 

/ \\+\\2\ (n+k-2i-2j-2l-l)/2 
_ (^)- t/2 c -H.-H/2. 



Also, 



(i + j-fc-Z)! 

Sfc+i-j [fc + Z-j]! / i- 1 



+ J -*-0 ! \k + l-j 

Finally, dominated convergence yields (7.2) and we are done. □ 

Theorem 4.1, for the case of the canonical process on S(M. ), will now 
follow immediately from Theorem 5.2 and the following result, modulo the 



GAUSSIAN KINEMATIC FORMULAE 



21 



issue of the moment assumption (5.7) which appears among the conditions 
of Lemma 5.2, but not among those of Theorem 5.2. We shall dispose of this 
issue in Section 8.2 below. 

Theorem 7.2. Suppose D satisfies the conditions of Lemma 7.1. Then 
in the notation of that lemma, and with JVQ(D) defined at (6.3), 

(7.4) p i {D) = {^T i Ml{D). 

Proof. We start by computing the p^. By Lemma 7.1, 

x Cj-i{p,ipkhi+j-k-lh-Di) 



(27T)-*/ 2 ]T (' ^M^P,^, 
m=Q ^ ' 



where the M are the generalized Minkowski curvature measures of (6.2). 
With the pi determined, we can now turn to the pj. By (5.9), these are 

L(i-1)/2J fiy 



L(t-1)/2J i-2l-\ 

'i-2l-\ 
rn 



M m+ i(D,tp k h 



i-2l-l-m) 



on applying the definitions of Hermite polynomials and AiJ{D). □ 

8. Beyond the canonical processes. In this section, we shall describe two 
different routes to extend the proof of Theorem 4.1 from the special case of 
the canonical process on the sphere to the more general Gaussian processes 
described in the theorem. The first route will indicate that this special case 
actually quite trivially implies many others. The second route is much longer, 
but gives the full result. 



8.1. The canonical process and processes with finite expansions. To save 
on notation, we shall for the moment consider only real valued processes /. 
The extension to the vector valued case is trivial. 
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It is well known that any continuous Gaussian process on a parameter 
space M has a representation of the form 

(8-1) = 

fc=i 

in which the are independent N(0, 1) and the <pk an orthonormal ba- 
sis on the reproducing kernel Hilbert space of /. Suppose that this expan- 
sion has I < oo terms, and that / has constant unit variance. Define the 
mapping tp:M— >M. 1 by setting ip(t) = (<pi(t), . . . , <pi{t)), and note that since 
1 = K{f 2 (t)} = £}=i f actually maps M to a set <p(M) C S(R l ). As- 

sume that <^9 is one to one, and define a new process / on ip(M), by setting 
for x £ <p(M), f(x) = f{p~ l {x)). Then 

l 

(8.2) E{f(x)f(y)} = J] = 

k=l 

That is, / is a version of the canonical Gaussian process on S(W). 

Now note that the excursion sets of / and / are related by the fact that 
A(f,ip(M),D) = (p(A(f,M,D)) and so, as long as ip is smooth enough (C 2 
suffices) it is easy to relate the Lipschitz-Killing curvatures of A(f, M, D) to 
those of A(f, ip(M), D). Their Euler characteristics, Co, for example, will be 
identical. As far as the others are concerned, the Lipschitz-Killing curvatures 
of a set A £ y(M), computed with respect to the usual Euclidean metric on 
S(M. 1 ), will be identical to those of ip~ l (A) € M computed with respect to 
the Riemannian metric on M which is the pull-back to M of the Euclidean 
metric on (p(M) by Another calculation like (8.2) shows that this is 
precisely the metric induced on M by the process /. 

These observations, along with the fact that Theorem 4.1 holds for the 
canonical process on S(W), now suffice to establish it for all processes satis- 
fying the assumptions of the theorem and for which the expansion (8.1) can 
be taken to be finite. 

8.2. Completing the proof for the general case. There is a general tech- 
nique in geometry that in the context of our problem, argues as follows: 
Suppose that we want to prove Theorem 4.1, in particular, (4.2). Suppose 
we could show that there exist functions p(i,j,D) dependent on i, j and D, 
but neither on the distribution of y nor on the topology of M, such that 

dim M—i 

(8.3) E{£ l (Mny- 1 (D))}= £ C l+j (M)p(i,j,D). 

3=0 

Then in order to identify the functions p(i,j,D), we could choose a param- 
eter space and random process that were simple enough to enable us to 
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compute E{£j(M r\y~ 1 (D))} in full. Writing it in the form of (8.3) would 
then allow us to determine the function p, and so we would have the result 
in full generality. In fact, this is precisely what we have done by working 
with the canonical Gaussian process on the sphere. 

An argument of this form is, for example, generally used to prove the 
kinematic fundamental formulae (1.6) and (5.2). 

All that remains, therefore, for us to have a full proof of Theorem 4.1 
is to prove (8.3). However, as is the case for the KFF, such representation 
theorems are difficult to prove, requiring considerable technicalities which 
go far beyond the level at which we have so far worked in this paper. The 
basic differential geometric tool that appears, and replaces the KFF that we 
have relied on so far, is the critical point theory of Morse. It is this, along 
with the smoothness assumption (4.1), that allows us to complete all the 
arguments. 

By agreement with the editor, who wanted to save space, and mindful of 
the fact that we would like to sell a few more copies of [3], we send you to 
Chapter 15 there, (which, in turn, will send you back to Chapter 13, which 
will demand that you also read Chapters 8-12) to see how this argument 
plays out in detail. 

Many other details and additional results that we left out of this paper 
for space reasons can also be found there. 
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